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Abstract 



We show that Haar-based series representation of the critical Riemann- 
Liouville process R a with a = 3/2 is rearrangement non-optimal in the 
sense of convergence rate in C[0, 1]. 
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Introduction and main result 

The aim of this note is to solve a problem stated by AAyache and W. Linde 
in their recent work [lj. In their article, the quality of several common series 
representations of fractional processes is considered. In particular, Ayache and 
Linde examine the representations of Riemann-Liouville processes R a , a > 1/2, 
based on Haar and trigonometric systems. They were able to show that Haar-based 
representation of R a is optimal w.r.t. the uniform norm when 1/2 < a < 3/2 while 
for a > 3/2 it is not optimal. Their approach however does not yield an answer in 
the delicate critical case a = 3/2. We will show that Haar-based representation is 
not optimal for R 3j/2 either. 

It is an immense pleasure for the author to stress that the main tool used in 
this note is due to V.N.Sudakov whose anniversary we celebrate in this volume. 

The reader can find recent results on series representations of fractional pro- 
cesses and fields in [2], [3], pE], [5], [6], [9], [TO]. 



*The work supported by the RFBR-DFG grant 09-01-91331 "Geometry and asymptotics of 
random structures" . 
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Recall that Riemann-Liouville process R a is defined by white noise represen- 
tation 



I r 

R a (t) = — — / (t - uY^dWiu), < t < I. 
Jo 

This process (RL for short) is known to have continuous sample paths whenever 
a > 1/2. It belongs to the family of so called fractional processes along with 
more widely known fractional Brownian motion W H (whenever a G (1/2, 3/2) and 
H = a — 1/2) and periodic stationary Weil process I a . All these processes differ 
by very smooth terms and therefore their approximation properties we discuss 
here are essentially the same, see [TJ for details. We also refer to [8] for further 
properties and stable extensions of RL-process. 

The representation of RL-process with a > 1/2 generated by Haar base writes 
as follows, 

oo ai-i 

Ra{t) = u f(^TT) + £ E (!) 

^ ' j=0 k=0 

where R a is the classical RL integration operator, 



1 /"* 

i2 Q /i(t) = — — / (t - u) a - l h(u)du, 

r(«) 7 



(Cj,fe)} is a family of i.i.d. standard normal random variables and hj^'s are 
the Haar functions 

hj k(t) — 2 J//2 \ lr 2k 2k+i \(t) — lrafc+i 2fe+2\(t) ?• . 

I [iJ+T'iJ+rJ [il+T'iJ+Tj J 

The series converges almost surely uniformly on t G [0, 1], i.e. it converges in the 
sense of the uniform norm 1 1 ■ W^. 

Recall that one can evaluate the integrated Haar functions by the formula 

, , w , I 312 [( 2k + 2\ a ( 2/fc + l\ a / 2k 



T(a + 1)[\ 2J+ 1 J + V 2i+1 )+ V 2i+1 

(2) 

Now we recall the necessary standard notation as well as the notions related 
to finite rank approximation of Gaussian random functions. 

Throughout the article f n >z g n means liminf n f n /g n > while /„ x g n means 
that both f n y g n and g n >z fn hold. We write #(-B) for the number of points in a 
set B, while \T\ denotes the length of an interval T. Finally, c denotes unspecified 
positive and finite constants which can be different in each occurrence. 
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Let X be a centered Gaussian random element of a normed space ■ ||). 

The ^-numbers (stochastic approximation numbers) of X are defined as 

n-1 "| 

Here the infimum is taken over all families of random variables and all finite 
deterministic subsets (<pi) C X. A series representation 

oo 

1=1 

is called optimal if 

x £ n (X) as n — > oo. 

It is called rearrangement non- optimal, if it can not be rendered optimal by any 
permutation of (<fi). In particular, for non-optimal representation the approxima- 
tion error E |E°^ n £« tends to zero slower than optimal rate l n (X). 

In the particular case of RL-process and uniform norm || • H^, the optimal 
approximation rate is well known, see [5J. Namely, for any a > 1/2 it is true that 

£ n (R a ) x n' {a - 1/2) Vh^l, as n -> oo. 

The optimal rate can be attained, for example, by using modified Daubechies 
wavelet base, see [10]. The remaining question is to understand which representa- 
tions provide this rate and which ones perform more poorly. 
For the critical case, a = 3/2, we have 

£ n (R 3/2 ) x n' 1 Vhin, as n -> oo. (3) 

Our main result shows that Haar-based representation is rearrangement non- 
optimal. 

Theorem 1 Let be any rearrangement of the Haar-based series represen- 

tation (JI]). Then 

y n~ Inn, as n — > oo. 



We observe a gap of order yum w.r.t. the optimal rate (J3j). 



t n {X) 



inf 



E 



E 



oo 



E 



5> 
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Proofs 

From now on, we consider only the critical case a = 3/2. 

Introduce the function H(t) = (t — 2)+ 2 — 2(t — 1)+ +t+ 2 which is the unsealed 
version of integrated Haar functions (j2J). Namely, we have 

(E h )(t) H ( 2J+1 (t-^r)) m 
Notice that for t > 2 

H{t) = ^J* (x 1 ' 2 -{x- l) 1/2 )dx ~ ^ t" 1 / 2 , as t -> oo. (5) 

For any level number j and a set of positive integers K C [0, . . . , 2 J — 1] we consider 
a Gaussian process {X K (t),t G [0, 2 J+1 ]} defined by 

X x (t) = ^ fc #(t-2A;) = &H(t-2k), 

keK keK,k<t/2 

where is a family of i.i.d. standard normal random variables. Clearly, Xk 
is the unsealed version of the i^-part of level j in Haar-based representation (pQ), 
since by (T4|) 

J2^ k ( i? 3/2%fc)(0 = 



23/2+ir(5/2) 



(6) 



in distribution. 



Proposition 2 There exists a numerical constant c such that 

E sup |Xx(t)| > c j, 

for all j > and all K C {0, . . . , 2-? - 1} such that #(K) > 2 j -\ 
Corollary 3 By scaling (0) it follows that 



E sup 

te[o,i] 

where hj^ are Haar functions. 



keK 



>cj2~i, (7) 
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Proof of Proposition [21 First of all, let us evaluate the natural distance 
associated to the process X K . Let < s < t < 2 J ' +1 . We have 

X K (t) = J2 &H(t-2k) 

keK:k<t/2 

£ k H(t-2k)+ Yl &H(t-2k) 

k£K:k<s/2 keK:s/2<k<t/2 

:= Y(s,t) + Z(s,t) . 

Notice that the variable Xk{s) belongs to span^, k e K, k < s/2) and the same 
is true for Y(s, t), while Z(s, t) is orthogonal to that span. By using (jSj), it follows 
that for s < t we have 

E{X K {t) -X K {s)) 2 > EZ{s,t) 2 = H{t-2k) 2 

k£K:s/2<k<t/2 

> c it -2k)- 1 . 

k€K:s/2<k<t/2-l 

Suppose now that we found in [0,2 J+1 ] an ordered subsequence t%, . . . ,t m , m > 
2i/ 2 - 2 such that for any 2 < i < m it is true that 

Y ft - 2k y l ^ c i- ( 8 ) 

fceif:ti_i/2<fe<ti/2-l 

Then, by Sudakov lower bound (see e.g. [7], Section 14) 



E sup \X K {t)\ > E sup X K (U 

te[0,2J+ 1 ] l<i<m 



21 1/2 



> c(lnm) 1 / 2 - inf [E(X(t h ) - X{t t2 )f] 

> cj 1/2 -cj 1/2 = cj. 

It remains to find a sequence (U) satisfying assumption (jSJ). Would K coincide with 

the maximal set, K — {0, . . . 2 J — 1}, we could just take tj = i2^ 2 , i = 0, . . . , 2 3 '/ 2 — 1. 
However, the situation is more delicate for general case. We will use the following 
elementary fact worth to be stated separately. 

Lemma 4 Let the points Si < s 2 < ■ ■ ■ < S2 q belong to an interval T in R. Then 

r-l 



«_ - «, - T H 



max _ 

l<r<2g ^—f S r — S{ \T\ 

i=l 
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Proof of Lemma 31 Take an integer j such that 1 < j < q. We have 

2q 2q r -l 

Y ( Sr - Sr -^ = Y Y ( si+i - 

r=q+l r=q+l l=r—j 

2q-l 

< Y j{si+i-s{)<j\T\. 

l=q+l-j 



By convexity of the function x — > - we have 



2q / 2q 

- Y - J — ^ 1 E 



-1 



[S r - S r _j) I > 



* r =q+l r r 3 X 1 r=q+l / v ! 

By summing up over j we get 

2q q q , 2q 

- Y Y—^~ = Y-Y 

$ r=q+l j=l Sr Sr ~i j=l q r=q+l S ' r Sr ~3 



> 



3=1 



Clearly, 

r—1 ) — 1 



max 



> max = max 



Kr<2q* — < S r — Sj q+l<r<2q * — ' S r — Sj q+l<r<2<j — ' S„ — S„_„- 

1=1 i=r—q j = l J 

1 2q 9 1 

q r=q+lj=l b r r -J 11 

We continue the proof of Proposition [2j For the sake of expression simplicity, 
assume that j is an even number, thus 2 jf//2 is an integer. We split the interval 
[0, V - 1] in 2- J '/ 2 blocks B L = [l2^ 2 , (l + 1)2^ 2 ) . 

Let 



J=| i: #(jfnfl,)>^/ , J 
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and m = We show that m is large enough, since 

2 J/2_1 

y- 1 < #(jo = E n 

t=0 

= J3#(ifns t ) + X;#(^nB t 



yields 



< mX' 2 + (2^ 2 - m) i 2^ 2 



lml 1 m 



thus 



2 - 4 4 2^/2' 

m > I 2 ^ 2 > 2^ 2 - 2 . 
~ 3 

For each block S t , t G I, we apply Lemma H] with T = B L , \T\ = 2^ 2 and 
{si, . . . , s 2(J } C P) B L and choosing g as large as possible, 

2q>#(KDB L )-l>^2i/ 2 -l. 

By Lemma HI we can find an integer point v L G {si, . . . , S2 9 } C X f] B L such that 



E — >% in( i> °j- 

s6iCn-B t ,s<v,, 



By renumbering the points {2v L , i G /} we get a large collection of m > 2 J / 2 2 
points ti with the property 



2 

k£K:ti-i/2<k<ti/2-l keK:ti-i/2<k<ti/2-l 



E fe-2*)- 1 = ^ E (^/2-^)- 1 

fcg_ft::t i _i/2<A:<t 4 /2-l 

> ^ E (t,/2 - A;)- 1 - 1 

k£K:ti-i/2<k<ti/2 

> \ E {v L -k)- l -l>c 3 , 



2 

keKf]B,,k<v L 

where the index i is such that tj/2 = v L . We also used here the fact that all t> t 's 
and all fc's are integers, thus we don't loose more than one while dropping one k 
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from the sum. Therefore, the collection {tj,z < m} satisfies ([S]) and Proposition [2] 
is proved. □ 

Proof of Theorem [TJ, Actually our theorem follows from Proposition [2] and 
Lemma 2.3 in [1] immediately. However, we recall the proof for reader's conve- 
nience. 

Let {£i(fi} be any rearrangement of the Haar-based series representation (pQ). 
Take an integer j > 1, let n = 2 J_1 and let K C [0, . . . , 2 J — 1] be the set of all k 
such that R a hj,k belongs to the set {tpi,i > n}. Clearly, =ff(K) > 2 J ' -1 , and we can 
write 

i>n keK 

where Y and the sum over K are independent. By standard arguments based on 
Anderson inequality for Gaussian processes (see [7]), we obtain 



E 



%>n 



> E 



k£K 



Now ([ZD yields 



and we are done. □ 



E 



%>n 



> c j 2 J x n 1 Inn, 



The author is grateful to A.Ayache and to W.Linde for interesting discussions 
and for attracting his attention to the problem solved in this note. 
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